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NGC 6441 is one of the most massive and most metal-rich globular clusters 
in the galaxy, and is noted for an unusual extended horizontal branch which 
reaches past the instability strip. We find evidence that there are two different 
populations of stars within the heavily-populated red clump. Once a differential 
reddening correction is applied, a large but compact group of stars is found at the 
faint red end of the clump in the color magnitude diagram. Brighter, bluer stars 
j^ . in the clump are found to be more centrally concentrated within the cluster at 

very high level of significance. Curiously, the blue horizontal branch stars show 
a more complex distribution, and are not more centrally concentrated than the 
O . brighter red clump stars. The spatial distributions of clump stars is in agreement 
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with the idea that the the brighter bluer part of the clump is a helium-enriched 
second generation. The blue horizontal branch stars may be showing evidence 



"^ . that they are being dynamically evaporated. 



Subject headings: globular clusters: individual: NGC 6441 — stars: abundances 
— stars: horizontal branch — stars: kinematics and dynamics 



1. Introduction 

Globular clusters (GCs) have long been seen as examples of simple stellar populations, 
though the current puzzle of multiple populatio ns began wi th the discovery of "chemical 



anomalies" within GCs over thirty years ago (see lKraftlll994j for a review of early studies). 
Many studies have shown that while generally homogeneous in Fe-peak elements, GC stars 
ca n exhibit internal sp reads in light element abundances (a more recent review can be found 
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Gratton et al.l 12004). These anomalies have been observed at both the turnoff and along 



the subgiant branch (JGratton et al.ll200ll ). and so the stars must have formed with at least 



some of the chemical anomalies. 
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High-precision photo metric data have pr ovided the means to identify broader evidence 
for multiple populations. iBedin et al.l ( 120041 ) used Hubble Space Telescope (HST) data to 
show that the main sequence (MS) in u Centauri is split into two components. S urprisingly, 



the bluer MS has twice the metallicity of the red component (jPiotto et al.l 12005), which can 



only b e explained by helium enrichment of the more metal- rich component. iPiotto et al. 
(120071 ) found that the MS of NGC 2808 is split into three components, and concluded this 
is probably also evidence of stellar populations with differing helium abundance. Splits on 
the subgiant branch (SGB) of multiple cl usters have a lso been reported, including u Cen, 
NGC 1851, NGC 2808, and NGC 6388 JPiottol boooj ). However, the oldest photometric 
evidence involves the unusual distribution of stars along the horizontal branch (HE). In this 
area, the clusters NGC 6388 and NGC 6441 are particularly noteworthy. Both clu sters have 
high metallicity (([Fe/H]) = —0.45 and —0.44, respectively; ICarretta et al.ll2009bl ). yet have 
extended horizontal branches that continue past the instability strip to T^ff > 30, ( 300 K . 
The RR Lyrae variables in both clusters have unusually long periods (jPritzl et al.ll2000l ). 
putting them in a third Oos terhoff class. These features may all be produced by populations 
of differing He enrichment JCaloi fc D'Antonal 120071 : IP'Antona fc Caloill2008h . 



While the HB of NGC 6441 has been studied several times before (lYoon et al.l 12008 



Busso et al.ll2007l : ICaloi fc D'AntonJ2007l : iRaimondo et ahlbooj iRich et aLlEoOTJ ) , two sub- 
jects have not been addressed: 1) although the signature of differential reddening in the 
cluster has been seen, researchers have not attempted to remove its effect on the CMD, and 
2) the spatial distribution of the stars throughout the cluster has not been examined. If a 
second generation of stars forms from gas enriched by a previous generation, and the gas 
has had a chance to cool and collect in the gravitational potential of the cluster, the second 
generation of stars is expected to have a more concentrated distribution. This ef fect seems to 



have been observed in several clusters, as in the two MS pop ulations of uj Cen (JBellini et al. 



20091). the giant branc h of NGC 3201 (ICarretta et al.ll2010l ). and potentially many others 
(ICarretta et al.ll2009al ). In a massive and old cluster, however, dynamical effects like mass 
segregation may affect the stellar distribution. We therefore set out to examine NGC 6441 
more closely. 



Data 



NGC 6441 is a massive bulge cluster with Mi 



V 



9.63 (JMcLaughlin fc van der Marel 



20051 ) . We rereduced HST WFPC2 archival data taken as part of proposal IDs 5667, 6095, 
8251, and 8718, and ACS HRC data taken in proposal ID 9835. These datasets, which 
contain optical and some near-UV observations, were selected in order to obtain nearly 
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complete samples of HB stars over a large fraction of the cluster's mass. The data cover 
from the center of the cluster to a ra dius of ~120", representing ~3 half-light radii [vh = 34'.'2; 
McLaughlin &: van der Marell l2005l ) . or about 15 core radii (r^ = 7'.'77). 



The WFPC2 images were analyzed using the HSTPhot package (JDolphiru 120001 ) . and 
stars were cross-identified among the multiple pointings. Median internal magnitude error 
estimates were 0.02-0.03 mag at the level of the horizontal branch for the images (from 
proposal 8251) covering the largest portion of the field discussed here. The ACS HRC 
images were processed using the ACS module of DOLPHOT (?). Populations were primarily 
identified via the color-magnitude diagrams (CMDs) shown in Fig. [1] — primarily ACS 
F555W versus F555W-F814W for the cluster core, and WFPC2 F439W versus F439W- 
F555W elsewhere. A different procedure was needed for fainter blue HB stars (as discussed 
below). The HB star group contains 2125, while the combined red giant/asymptotic giant 
(RGB/AGB) group 912 stars. The HB can be subdivided into the red clump (RC; 1754 stars) 
and the "blue" HB (BHB; 371 stars). The identification of the BHB sample deserves some 
addi tional explanatiori . T he sample described here includes RR Lyrae variables (identified 
from lPritzl et al.ll2003l and lCorwin et al.ll2006l ) and some red HB stars bluer than the clump, 
but does not include extreme HB stars (771555 > 20) as those stars were found to have 
significant incompleteness. The BHB stars were validated using near-UV data from proposal 
8718 (PI Piotto) and 5667 (PI Margon; F336W filter only) when available because crowding 
of star images is reduced in the UV where cooler giants and main sequence stars emit less. 
UV data was only available for two mostly overlapping WFPC2 pointings, and so did not 
cover the entire extent of the optical data. Importantly though, it does cover the center of 
the cluster where blending with brighter cluster stars is the most likely to cause BHB stars 
to be lost in optical photometry. In the F255W filter for example, the blue HB stars are 
among the brightest cluster sources, and crowding effects are negligible. Outside of the fields 
covered by the UV images, the F439W filter plays a similar role by enhancing the brightess 
of BHB stars relative to the numerous faint RGB and MS stars. Our artificial star tests 
indicate that incompleteness among our BHB sample is less than 3% in the central regions. 



Nataf et al.l (1201 ll ) and lMilone et al.l ( 120121 ) have shown that structure within the RC of 



47 Tuc can also reflect population differences within a cluster. Because the RC population in 
NGC 6441 is large, it provides an excellent probe of the radial structure of the cluster even 
when subdivided. In light of this, we examined several different ways of bisecting the RC 
population: by magnitude, by line parallel to the slope of the HB, and by line perpendicular 
to the slope of the HB. In this way we can look for signs of population differences somewhat 
independent of theoretical expectations. Because the results of sub dividing the clump cou ld 
be affect ed by the cluster's differential reddening (see Fig. 7 of iRaimondo et al.l |2002| or 
Fig. 4 of IBusso et al.l 120071 ). we corrected for this by determining the median color shift of 
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lower RGB stars (18.3 < 7/1555 < 19.2) from a mean line in the CMD. The correction was 
determined for each HB star from RGB stars within 9" or 15" depending on the number of 
RGB stars available. In the center of the cluster, we relied on the high resolution photom- 
etry from the ACS HRC (proposal 9835, PI Drukier) to reduce incompleteness and also to 
minimize the effects of blending in the photometry. If uncorrected, this allows some clump 
stars to be measured as too bright and allows some red giants to masquerade as clump stars 
when blended with fainter stars. Our artificial star tests indicate that our incompleteness 
for clump stars is less than 1% across the entire i7i?C field, and the median magnitude offset 
between input and output stars was approximately 0.002 mag. 

Sky coordinates were determined using the METRIC task in the STSDAS IRAF package 
Before examining the radial distributions of star s in the cluster, we recalc ulated the cluster 
centroid using a method similar to that used by iMontegriffo et al.l ( 1995[) , and found to be 
a = 17^50"1?:9 and 6 = -37 : 03 : 05.93, or about 6" from that of lHarriJ Jl996h . 



3. Results 

We used Kolmogorov-Smirnov (KS) test statistics to calculate the probabilities that 
the cumulative radial distributions (CRDs) of different subpopulations were drawn from the 
same underlying distributions (Fig. |2]). In our analysis, the RGB/AGB group was used as 
a stand-in for the cluster as a whole as this portion of the CMD does not show multiple 
populations as clearly. 

When examining subdivisions of the red clump, we found that cuts parallel to the HB did 
not produce radial distributions of the bright and faint components that were significantly 
different. Cuts perpendicular to the HB and cuts by F555W magnitude both resulted in 
groups that had significant differences in their radial distributions, and the results of KS 
tests were essentially equivalent. (This is perhaps expected given the steep slope of the HB 
at the location of the RC.) For the lack of a strong reason for selection, we concentrate 
on cuts made at a F555W magnitude of 17.9 for the rest of this study, and designate the 
subpopulations the bright red clump (bRC) and faint red clump (fRC). As seen in Fig. |2l 
the bRC is quite significantly more centrally concentrated than the fRC, with the probability 
that the two are drawn from the same underlying distribution being 8 x 10^^^. The bRC is 
also more centrally concentrated than the RGB/AGB. 

A closer examination of the clump reveals clues to why this is. After the clump has 
been corrected for differential reddening, the color extent of the clump is similar in the three 
radial bins we selected, indicating that the issue has been largely removed. In the innermost 
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bin, there are few stars at the faint red end of the clump (^555 ^ 17.95). There is a strong 
concentration of stars at that end of the clump in the middle radius bin, but in the outermost 
bin the stars are distributed nea rly uniformly acros s the clump. This behavior can be seen 
in the CMD array in Fig. 4 of iBusso et al.l ( 120071 ) as well — there is a tendency for the 
clump to appear stubbier and more concentrated in color outside the centralmost part of the 
cluster. 

The BHB shows yet another pattern in its radial distribution. Overall, the CRD (Fig. [2]) 
indicates that the BHB is less centrally concentrated than both the bRC and the RGB/AGB. 
Comparison of the BHB and RGB/AGB distributions produces a KS probability (that they 
derive from the same underlying radial distribution) of 4 x 10~^. But it is also interesting 
that the slope of the distribution between ~25" and ~60", is lower than those for other 
populations. The shallower slope perhaps suggests that the radial distribution of BHB stars 



is bim odal, though more weakly than is often seen in blue straggler studies (JFerraro et al. 
20041 . e.g.). 



To clarify the situation further, we have compared the stellar populations with the clus- 
ter light distribution using do ubly-normahzed rati os Rpop = {,Npop/Nl°^f )/{L sampled / L^slmpied) 



from blue straggler studies ( iFerraro et al.l Il993[ ). We calculate d the sampled luminosity 
L,^„mpi e.d using a King model with structural parameters taken from lMcLaughlin fc van der Marel 
( I2OO5I ) that fits the surface brightness distribution of the cluster. The ratio is calculated in 
annular radial bins centered on the cluster, and is designed to equal 1 if the stellar popu- 
lation matched the light distribution (see Fig. [3]). In the Rpop ratios, the BHB stars show 
an enhancement in numbers in the third radial bin (60 — 90") that is significant at the 2a 
level. The bRC stars show the most significant deviations from following the light, with the 
population ratio steadily decreasing moving outward from the center. 

Before discussing the physical meaning of these distributions, it is worth commenting 
on possible complicating factors. Although others have concluded that differential reddening 



is not responsible for the slope of the the clump ( IBusso et al.l 120071 : iRaimondo et al.l 12002 



Sweigart &: Catelanlll998l ). we applied a correction as described in §2 to eliminate the ma- 
jority of its effect. The referee wondered whether blending could produce the strong central 
concentration seen in the bRC population. This idea is contradicted by the fact that strong 
central concentration among the brighter clump stars is still seen when the clump is split 
along a line perpendicular to the HB (more like a color cut) while the concentration disap- 
pears when the cut is made parallel to the HB (more like a magnitude cut). It could also be 
imagined that incompleteness might be an issue in the population of BHB stars (particularly 
among the faint blue tail in the optical CMD), but there are several ways to see this is a 
minor consideration. First, in the core where incompleteness is the biggest issue, our star 
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selections made use of near-UV photometry — this reveals the hottest HB stars even when 
very close to redder HB stars and giants, and crowding issues are negligible in the near UV. 
Second, the blue tail of the BHB is a small fraction (about 13%) of the stars in the BHB 
sample. So while we may have lost some (~3) of the hottest HB stars (faintest in the optical) 
to incompleteness, this would not affect conclusions about the radial distributions. 



4. Discussion 

The impression from the red clump stars in NGC 6441 is that the bluest and brightest 
of the red clump stars are considerably more centrally concentrated than the stars at the 
faint red end of the clump, which are found in a small area of the CMD. However, the 
even hotter BHB stars are less centrally concentrated than the bRC stars and show a radial 
distribution of different character. These indications of three HB populations with differing 
radial distributions in NGC 6441 are difficult to fit into a consistent picture. 

Cluster self-enrichment is the current leading explanation for multiple populations in 
globular clusters, involving a second generation of stars formed from gas polluted in he- 
lium and some heavier elemen ts by ejecta from a group o f stars such as intermediate-mass 



asymptotic giant bran ch stars fIVentura &: D'Antonall2008l ) or rapidly-rotating massive stars 



(iDecressin et al.l 120071 ) in the first generation (FG). There are several factors that affect the 
radial distribution of the second generation (SG), although one of the primary predictions 
is that the formation of the second generation only begins once gas has cooled and concen- 
trated in the potential well of the cluster. SG stars are thus expected to be more centrally 
concentrated than the FG, although later dynamical evolution of the cluster could modify 
this. 



Caloi fc D'Antonal (120071 ) produced synthetic models of NGC 6441 HB stars using this 
picture with the fRC representing FG stars (helium abundance Y = 0.25, and having a 
compact distribution in the CMD), while the bRC contained SG stars reaching helium abun- 
dances of at least Y = 0.35. BHB stars in the models involved a combination of stars with 
even higher abundances (up to y = 0.40) and stars that had evolved from the bright end of 
the clump (large helium abundances result in large blueward excursions during HB evolu- 
tion). They estimated that the enriched HB population contained 60% of the stars, although 
they based their study on a smaller sample of stars and did not account for differential red- 
dening, which smears out the peak at the faint end. This last factor also influences whether 
or not a gap in Y between FG and SG stars is inferred (see their Fig. 7) 

On the basis of the CMD, there is good evidence of two populations. However, the 
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radial distributions paint a more complex picture. Specifically, the BHB stars do not not 
continue the trend that hotter HB stars (presumably even more helium enriched) are also 
centrally concentrated. 

Though the second generation may be born with a more centrally concentrated distribu- 
tion, simple physics tells us two-body interactions must occur, causing the SG stars to diffuse 
outward over time. This effect would be enhanced by the fact that helium-enriched stars are 
expected to be less massive th an stars of normal compo sition because they can evolve off the 
main sequence more quickly (JD'Antona fc Caloil 120081 ) . Two-body interactions in a cluster 
lead the stars toward equiparti tion of kinetic energy, "heating" the SG stars and giving them 
a more extende d distribution (ICarretta et al.ll2009al ). This process is expected to act within 
a Hubble time ( iD'Ercole et al.ll2008[ ). affecting the innermost part of the distribution most 
quickly. Over time, the cluster core should begin to be depleted of the progenitors of BHB 
stars, while stars on larger orbits would be less affected. 

As an exercise, we can calculate the radius at which the dynamical relaxati on timescale 



equals the lifetime of the cluster. We ag ain use cluster structural pa rameters from lMcLaughlin fc van der M 
( I2OO5I ) and the following equation from iBinney fc Tremaind (J2008[ ) , which describes a stellar 
system that is a singular isothermal sphere: 



t 



relax 



5 Gyr IMq a 

In A m 10 km s"^ V^PC 



and find that the relaxation radius is ~60". This falls within the part of the BHB radial 
distribution that has shallower slope than the rest of the stellar groups. Qualitatively, the 
BHB distribution shows the signatures of dynamical modifications. 



Other mechanisms such as gas expulsion (JDecressin et al.ll2010l ) can affect the dynamical 
evolution of the cluster, so that more detailed models will be needed to verify whether the 
radial distributions are consistent with expected masses for helium-enriched populations. 
Examination of other dynamically interesting populations, such as blue stragglers, would also 
clarify the state of the cluster structure. In the more massive blue straggler star populations 
of many clusters, there is a bimodal radial distribution, where the "zone of avoidance" 
corresponds well to the area expected to be cleared out dynamical friction. However, in 
NGC 6388, a cluster similar to NGC 6441 in many resp ects, blue stragglers still populate 
that area, making the cluster appear dynamically young (JDalessandro et al.l 120081 ). 



We thank the anonymous referee for very helpful comments on the paper. This research 
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Fig. 1. — The CM diagrams along the top show the populations used in the analysis. These 
are the faint red clump (fRC; green), the bright red clump (bRC; light blue), the blue 
horizontal branch (BHB; dark blue) and red giant branch (RGB; orange). The leftmost 
CMD shows data from the ACS HRC while the other two show WFPC2 data. The three 
CMDs along the bottom are zoomed in on the red clump at three different radial bins. The 
innermost bin uses HRC data and the others WFPC2. 
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Fig. 2. — Cumulative radial distributions of the analyzed populations. 
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Fig. 3. — Radial distribution of the double-normalized ratios for the selected populations. 
Black crosses mark the boundaries of the chosen annuli. The points are offset horizontally 
from one another for clarity. 



